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ABSTRACT 

Chandra X-ray Observatory grating spectra of the supergiant X-ray Binary 
4U 1700-37 reveal emission lines from hydrogen and helium-like S, Si, Mg, and 
Ne in the 4-13 A range. The spectrum also shows fluorescent lines from S, Si, 
and a prominent Fe Ka line at 1.94A. The lines contribute to the previously 
unaccounted "soft excess" in the flux in this range at orbital <p ~ 0.7. The 
X-ray source was observed during intermittent flaring, and the strengths of the 
lines vary with the source state. The widths of the lines (FWHM « 1000 — 
2000 km/s) can result from either Compton scattering or Doppler shifts. Power 
spectra of the hard X-rays show red noise and the soft X-rays and lines show in 
addition quasiperiodic oscillations (QPOs) and a power-spectral break. Helium- 
like triplets of Si and Mg suggest that the gas is not in a pure photoionization 
equilibrium. We discuss whether resonant scattering could affect the line ratios 
or whether a portion of the wind may be heated to temperatures T ~ 10 6 K. 
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1. Introduction 

High Mass X-ray Binaries (HMXBs) consist of a compact object (a neutron star or 
black hole) accreting from an early type stellar companion. Stars of spectral types O and B 
typically drive stellar winds with mass loss rates of 10~ 7 — 10~ 5 M Q yr _1 , and a portion of 
this outflowing gas may be captured by the gravity of the compact object. Simple analytic 
expressions of gravitational capture (e.g. Bondi & Hoyle 1944) show that in some systems, 
the observed X-ray luminosity may be accounted for entirely by accretion of the wind. Other 
HMXBs, however, are probably powered by Roche lobe overflow leading to a gas stream and 
accretion disk. Even systems in which it is thought that wind accretion dominates are 
often close to filling their Roche lobes (Petterson 1978). As the surface of the primary star 
approaches the critical potential surface, there is probably a smooth transition between a 
stellar wind enhanced along the line between the stars and Roche lobe overflow (Friend & 
Castor 1982). 

One of the earliest tools for understanding these systems, suitable for low resolution 
X-ray spectroscopy, was the examination of the orbital variation of the X-ray absorption. 
Buff & McCray (1974) pointed out that as an accreting compact object orbits a normal 
companion with a spherically symmetric stellar wind, the wind should absorb soft X-rays at 
all phases of the orbit. The amount of absorption should increase as the compact object falls 
behind the densest parts of the wind, close to eclipse ingress and egress. This technique has 
been used, for example, by Haberl, White, & Kallman (1989, hereafter HWK), who observed 
a binary orbit of 4U 1700-37 using EXOSAT. In this system, however, the absorption is not 
symmetric about orbital phase <fi = 0.5 (Branduari et al. 1978). HWK showed that the 
excess absorption at <fi > 0.6 could be explained by an accretion stream conveying 8% of the 
gas that accretes onto the compact object. 

In addition to a gas stream, an "accretion wake", a bow shock trailing the compact 
object, has also been invoked to explain X-ray absorption in several systems (Livio, Shara, 
& Shaviv 1979) and has been simulated numerically (Blondin et al. 1990). 

Another structure the compact object may create within the wind is a "photoionization 
wake", which results when the accelerating stellar wind crashes into gas ionized by the X- 
ray source. The ionized gas coasts, as it lacks the ions that accelerate the wind through 
resonance scattering of the stellar continuum (Fransson & Fabian 1980). Photoionization 
wakes are probably responsible for optical absorption lines from the primaries in Vela X-l 
and 4U 1700-37 at > 0.5. A gas stream or accretion wake could cover the compact object 
but not enough of the primary to explain the optical absorption (Kaper, Hammerschlag- 
Hensberge, & Zuiderwijk, 1994). Photoionization wakes appear in numerical simulations of 
X-ray sources embedded in stellar winds for a wide range of X-ray luminosities (Blondin 
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1994). 

High resolution X-ray spectroscopy is emerging as a powerful diagnostic of conditions 
in HMXBs. An ASCA eclipse spectrum of the X-ray pulsar Vela X-l showed emission lines 
that could be attributed to hydrogen-like and helium-like ionization stages of heavy elements 
(Si, Mg, and Ne, and possibly S, Ar, and Fe) as well as neutral fluorescence lines. The line 
strengths could be qualitatively understood through a model of photoionization of the stellar 
wind (Sako et al. 1999). ASCA observations of the HMXB pulsar GX 301-2 also showed 
fluorescent lines from Ne, Si, S, Ar, and Ca, as well as a 0.8 keV bremsstrahlung component 
that may be caused by shocks trailing the neutron star (Saraswat et al. 1996). 

The new generation of X-ray telescopes, Chandra-AXAF and Newton-XMM, are con- 
firming and expanding on the results suggested by ASCA and are revealing diverse new 
phenomena. Paerels et al. (2000) observed emission lines and radiative recombination fea- 
tures in the HMXB Cygnus X-3 with Chandra. Line ratio diagnostics suggest the emission 
lines are formed by X-ray photoionization of the stellar wind. Schulz et al. (2002b) observed 
the Vela X-l lines with much higher resolution than with ASCA using Chandra, and resolved 
lines that could be used to diagnose the state of the emitting gas. The HMXB black hole 
candidate Cygnus X-l shows variable emission and absorption features (Schulz et al. 2002a), 
and may reveal the presence of a wind focused by the gravity of the compact object (Miller et 
al. 2002). The superluminal jet source GRS 1915+105 shows absorption from H and He-like 
Fe (Lee et al. 2002), while SS 433 shows relativistically Doppler-shifted emission lines that 
provide density and temperature diagnostics for the gas in the jets (Marshall, Canizares, & 
Schulz, 2002). 

A new dimension can be added to high resolution X-ray spectroscopy of HMXBs when 
line features are observed through times when the X-ray continuum is varying. Power spectra 
of X-rays from HMXBs typically reveal noise with a fractional rms variability of 20-30% 
(Belloni & Hasinger 1990). Quasiperiodic oscillations (QPOs) with mHz frequencies have 
been observed from X-ray pulsars such as LMC X-4 (Moon & Eikenberry 2001a) and Hercules 
X-l (Moon & Eikenberry 2001b, Boroson et al. 2000). Observations of X-ray lines from a 
stellar wind with variable X-ray illumination can test photoionization models. 

For this investigation, we focus on the supergiant X-ray binary 4U 1700-37, discovered 
with Uhuru (Jones et al. 1973). The X-ray source is eclipsed with a 3.41 day period by 
its companion HD 153919 (Penny et al. 1973, Hutchings et al. 1973), determined to be an 
06.5Iaf star of M Q = 52 ±2 M and T cS = 42, 000 K. Doppler shifts of lines from HD 153919 
indicate an orbit that is nearly circular (Heap & Corcoran 1992). EXOSAT and BATSE have 
set an upper limit of ~ 4% on coherent pulsations near the expected ~ 100 s period of a 
typical X-ray pulsar (Doll & Brinkmann 1987, Rubin et al. 1996). 
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Corbet (1986) showed that most HMXB X-ray pulsars are either orbiting Be stars, or 
fall into separate period classes, 10 2 — 10 3 s and 1-10 s, that may represent wind-fed and 
disk-fed systems, respectively. The few reported detections of pulsations in 4U 1700-37 have 
not been repeated. Murakami et al. (1984) reported 67.4 second pulsations during a flare 
observed with Tenma. 

There is no consensus on the cause of the flares in 4U 1700-37. The X-ray light curve 
shows periods, lasting on the order of an hour, when the flux increases by up to a factor of 
~ 10. These flaring periods are accompanied by flickering in the light curve, also by about 
a factor of 10, on time scales of seconds to about 10 minutes. Brinkman (1981), based on a 
model of an entirely wind-accreting system, suggested that the flares were associated with 
accretion from the magnetotail of a neutron star. Tentative detection of pulsations suggested 
to Muakami et al. (1984) a similarity between 4U 1700-37 and LMC X-4, an X-ray pulsar 
that increases the fraction of pulsed X-rays during flares in which the X-ray luminosity can 
reach several times the Eddington luminosity. Levine et al. (2000) suggested three possible 
causes for the flares in LMC X-4: an accretion instability, fusion on the neutron star, or 
magnetic reconnection. In contrast to LMC X-4, in which the spectrum becomes softer 
during flares, the spectral changes during flares in 4U 1700-37 are not consistent, although 
there is evidence that the soft X-ray excess can increase during flares (HWK). 

The lack of a clear detection of pulsations in 4U 1700-37 is surprising, as most neutron 
stars in HMXBs pulsate. The neutron stars in LMXBs are rarely pulsars, but they may 
have weaker magnetic fields owing to physical effects of accretion over the longer life of the 
system (Bhattacharya & Srinivasan 1995). 

Measurements of the mass of the compact object in this system have been contradictory. 
The mass can be determined from the mass function of the O star, 

f{M ) = PK*(1 - e 2 f 2 /2nG = M x sm 3 i/(1 + l/q) 2 (l - e 2 f' 2 (1) 

where K Q is the orbital velocity of the O star, P is the orbital period, i is the inclination, 
q = M x /M is the mass ratio, and e is the eccentricity of the orbit. For K in km s -1 and 
P in days this gives: 

f{M ) = 1.038 x l(T 7 iir&P(l - e 2 ) 3/2 ( 2 ) 

Heap & Corcoran (1992) using Ko = 18 ± 3 km s _1 , as determined from IUE obser- 
vations, find f(M ) = 0.0023 ± 0.0005 M , resulting in an M x = 1.8 ± 0.4 M . Stickland 
& Lloyd (1993) re-examined the IUE spectra that led to this determination, recalibrated 
the wavelength scale, and found instead K = 5 — 10 km s _1 , leading conservatively to 
f(M ) = 0.00035 M , providing an uncomfortably low mass for the compact object. Stick- 
land and Lloyd raise the possibility that the measured UV line variations do not give the 
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radial orbital motion of the O star. Recently, Clark et al. (2002) using a measurement of 
Kq = 20.6 ± 1.0 by Hammerschlag-Hensberge et al. (work not yet published that involves 
reanalysis of spectroscopic data) found M x = 2.44 ± 0.27 M , uncomfortably high for a 
neutron star but lower than expected for a black hole candidate. 

We list the measured parameters of the system in Table 1, including the range of mass 
determinations for the compact object is in a close orbit (asini = 2.0 ± 0AR o , where 
Ro = 18 ± 3R , Heap k Corcoran 1992). 

2. Observations 

The Chandra observations we report here took place on August 22-23, 2000. We sched- 
uled the observation to coincide with the period before eclipse at = 0.65 — 0.8, so that we 
could investigate the pre-eclipse absorption component and the soft excess. The Chandra 
observation was coordinated with simultaneous observations with the Rossi X-ray Timing 
Explorer (RXTE, Jahoda et al. 1996). 

We used the High Energy Transmission Grating Spectrometer (HETGS, Canizares et 
al. 2000) aboard the Chandra X-ray Observatory (Weisskopf et al. 2002). We used the 
telescope's "faint" data mode, with a timed readout every 3.24104 seconds. The HETGS 
includes both a Medium Energy Grating (MEG) and a High Energy Grating (HEG) giving 
overall a spectral range of 1-35A, with a resolution of A/AA = 1400 at 12 A and A/AA = 180 
at 1.8 A. The dispersed spectra were then detected with the ACIS-S linear array of CCD 
chips. The zeroth order undispersed spectrum was severely piled-up, and was not used in 
our analysis. 

The observation was uninterrupted with a total exposure of 4.241 x 10 4 s, including a 
deadtime correction of 1%. Including photons in both the MEG and HEG, but excluding 
background, we detected 630,787 dispersed photons. 

2.1. Light Curves and Spectra 

We show a light curve for the combined HEG and MEG first order spectra in Figure 1. 
We have ignored photons near the wavelength of the Fe K line, where pileup may flatten the 
light curve during flares. 

We have determined the orbital ephemeris from RXTE ASM data to be T = 50088. 625± 
0.030 MJD, P = 3.411625 ± 2 x 10~ 5 days, and P/P = -1.5 ± 0.5 x 10~ 6 yr" 1 . 
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Along with the Chandra light curve, we also plot in red the 2-9 keV light curve as 
observed with the Proportional Counter Array (PCA) on board RXTE. We show only counts 
from the PCU 3 unit, which is the most reliable, although three PCUs were operating during 
most of the observation. The agreement between the Chandra and RXTE light curves is 
generally excellent, although the RXTE light curve appears systematically brighter near 
flare peaks. The RXTE observation extended 41.4 hours before the Chandra observation 
and 6.8 hours after. The peak count rate in the RXTE PCU 3 over the entire RXTE 
observation was about twice the peak count rate during the Chandra observations. 

For further analysis we divided the Chandra observation period into three main segments 
(Figure 1), which we call Flaring (F), Quiescent (Q), and Ending (E). The Flaring period 
starts at the beginning of the Chandra observation and ends at MJD 51778.6296, when the 
Quiescent period begins. (These times are spacecraft times, not corrected to the solar system 
barycenter.) The Quiescent period ends at MJD 51778.7454, when the Ending period begins. 

The division of the light curve into F, Q, and E periods is made empirically based on the 
count rate. The average Chandra count rates in the (F,Q,E) periods are (6.4,0.7,1.7) cts s _1 . 
However, as we discuss in §3, the absorbing column increases throughout the observa- 
tion, so that the intrinsic luminosity is more similar during the F and the E periods than 
the count rate would indicate. From simple broad-band fits of the spectra to absorbed 
power law models (see §3.1), we find average X-ray luminosities in the (F,Q,E) periods are 
(1.1,0.2,0.8) xlO 36 erg s -1 (corrected for absorption). 

In Figure 2, we show the l-lOA Chandra spectra of 4U 1700-37 at each of the three 
periods of the observation. We have co-added HEG and MEG counts in the ±1 orders. 
The background is negligible and has not been subtracted. We have applied standard CIAO 
procedures for "destreaking" the ACIS chips, and have verified that the strong lines do not 
overlap with chip gaps. Although for display purposes we show the total counts in the ±1 
orders of the MEG and HEG added together, we fit the models to each order and instrument 
in parallel, using corresponding response files. Fluorescence lines, notably an Fe K line near 
1.94A, are prominent in both the F and E periods, and less prominent in the Q period. 

High ionization lines due to H-like and He-like stages are also prominent in the 4- 12 A 
wavelength range. We identify these lines using House (1969), the analysis of lines in the 
Vela X-l spectrum by Schulz et al. (2002b), and the XSTAR code (Bautista & Kallman 
2001, Kallman & Bautista 2001). 
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3. Analysis 

To process and select counts, we used CIAO ( Chandra Interactive Analysis of Observa- 
tions) software version 2.2.1. 

We fit the continuum using the ISIS software package (Houck &Denicola 2000) and re- 
binning the spectra by factors of 20 so that each bin has enough counts to allow use of the 
X 2 statistic. This analysis is described in §3.1. 

To fit the spectral lines, we used the Sherpa spectral fitting tool. Instead of x 2 5 we use 
the C statistic (Cash 1979), which is appropriate when bins have few counts. We analyze 
the fluorescence lines in §3.2 and the high ionization lines in §3.3. 

For the spectral response files, we used acisheglD1999-07-22rmf N0004 . f its for the 
HEG and acismeglD1999-07-22rmf N0004 .fits for the MEG. We created auxiliary response 
files (ARFs) using the CIAO task "fullgarf" for event files for each of the three periods of 
the observation (Flaring, Quiescent, and Ending). 

3.1. Continuum Spectrum 

The hard continuum was fit using a power law absorbed by cold gas (a "Wisconsin 
absorber", Morrison & McCammon 1983). We need a second component to fit the soft 
X-ray continuum. We obtained adequate fits when we allowed for a second power law 
component with the same photon power law index but with a different normalization and 
behind a different column density. Such a component could arise from Compton scattering 
in the wind. Haberl & Day (1992) fit Ginga spectra of 4U 1700-37 to a similar model, and 
Haberl (1991) performed Monte-Carlo simulations of scattering of the 4U 1700-37 spectrum 
in a stellar wind. 

Parameters describing the continuum fits to the F, FH, FL, Q, E, and EL states are 
shown in Table 2. 

The average spectrum during the Ending period with our determination of continuum 
level and identification of spectral lines is shown in Figure 4. 

3.2. Fluorescent Lines 

We discuss the low-ionization fluorescence lines (Fe K, Ca, Ar, S, Si, and Mg, in order 
of increasing wavelength) observed with Chandra, paying special attention to the brightest 
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fluorescence feature, the Fe K line near 1.94 A. In Figure 3 we show the spectrum in the 
region of the Fe K lines during the F, Q, and E periods, along with our model fit. The fit 
in Figure 3a of the Fe K lines during the Flaring period shows that the data are slightly in 
excess of the model. We note that the model is not fit to the summed data, but that the 
spectrum in each HETG arm is fit in parallel. The HEG +1 and —1 spectra show Fe K 
peaks offset by 0.005A from each other 

The Fe K line may be subject to some pileup during the Flaring period. Two pho- 
tons with energies E\ and E 2 striking the detector within a frame time register with an 
energy E\ + E%. While the number of piled up photons can in some cases be deduced from 
the ACIS-detected energy, the energy carried by piled-up Fe K photons would exceed the 
ACIS threshold. Pileup in the Fe K line in the Flaring period would not only lead to an 
underestimate of the flux in the line, but also an overestimate of the line width. 

The edge apparent near 1.74A in Figure 3c results from the neutral Fe K edge introduced 
by the "Wisconsin" model of interstellar absorption we use to fit the continuum (wabs). 

Fitting the Flaring period spectrum to an absorbed power law plus Gaussian emission 
line gives a line centered at 1.937±0.002A (all errors are 90% confidence and are entirely 
determined from fitting) with a Full Width at Half Maximum (FWHM) of 2200 ± 700km s _1 . 
During the Quiescent period, the centroid of the line is at a marginally longer wavelength, 
1.940±0.002A, with only an upper limit of 2200 km s 1 for the line width. During the 
Ending period, the line is centered on a wavelength of 1.936±0.00lA but with a significantly 
narrower FWHM of 900 ± 400 km s -1 . 

To identify the ionization stages of Fe responsible for the line, we use the calculated Ka 
line wavelengths of House (1969). We find that the 90% confidence range is compatible with 
Fe II XII. 

The most likely location for emission of the Fe K line is the stellar wind. HWK point out 
that the strength of the Fe K line seen in mid-eclipse with EXOSAT is consistent with the 
out of eclipse flux and the portion of the wind blocked by the normal star or not illuminated 
by X-rays during eclipse. However, the possibility remains that some fluorescent emission 
arises on an accretion disk or gas stream. 

Motion of the stellar wind can introduce wavelength changes via the Doppler shift. 
Based on high resolution UV spectra obtained with IUE, Heap & Corcoran estimate a wind 
terminal velocity of = 2100 km s _1 , while van Loon et al. (2001) estimate a terminal 
velocity of Voo = 1700 km s -1 and a fractional microturbulent velocity of o = 0.15, which 
leads to a similar Doppler broadening. However, we are not able to separate the effects of 
ionization stage, Doppler shift, and Compton scattering in the appearance of the emission 



- 9- 



lines. 

Comptonization can affect both the apparent wavelength of an emission line and its 
width (Kallman & White 1989). Ignoring the weak temperature dependence in the Comp- 
tonization broadening, the equations of Kallman & White, based on models of Podnyakov, 
Sobol, & Sunyaev (1979), Ross, Weaver, & McCray (1978), Langler, Ross, & McCray 
(1978), and Illarionov et al. (1978), imply that AA/A = 0.0292r Th where r T h is the op- 
tical depth for Thompson scattering. This implies r^h = 0.25 ± 0.08 for the flaring period 
and Txh = 0.10 ±0.05 for the ending period. These values should be considered upper limits, 
as some line broadening could also result from a range of ionization stages in the emitting 
Fe or from Doppler shifts in the stellar wind, if the Fe K line is emitted in the wind. 

Assuming such a Txh, and that the line centroid wavelength A is shifted by Comptoniza- 
tion from the rest wavelength Ao according to X/X = 1 + 0.0375rxh (as given by Kallman & 
White, ignoring the temperature correction), we find that our observations are compatible 
with all stages of Fe up to Fe XVIII. We emphasize that these limits are set by the centroid 
of the Fe K line and by its observed width. 

The greatest scattering optical depth found from the line width, rxh = 0.25±0.08, found 
during the Flaring period, corresponds to a column density Nh = 3.8 ± 1.2 x 10 23 cm -2 . 
The absorbing column density seen throughout the orbit varies from « 4 x 10 22 cm -2 from 
= 0.2 to (p = 0.6, and 2 - 4 x 10 23 cm" 2 from = 0.8 - 0.9 (HWK 1989). Although 
EXOSAT did not observe from = 0.6 — 0.75, the models of HWK predict a range of 
column densities of N H = 0.4 — 3 x 10 23 cm -2 . The scattering optical depth r T h decreases to 
< 1.5 ± 0.8 x 10 23 cm~ 2 in the Ending period. Between the Flaring and Ending periods, the 
Fe K line is scattered by a decreasing column density, in anticorrelation with the absorbing 
column density, which increases between Flaring and Ending periods. This could be the 
result of the particular scattering geometry involved. The absorbing column density samples 
the wind along the line of sight from the X-ray continuum source, whereas the Fe K line is 
emitted in a more extended region in the wind (it persists through eclipse.) 

During the F period, fits show only a marginal Fe K/3 line. However, the statistics are 
better during the E period, when the lines are narrower. We fit the K/3 line to a gaussian 
with central wavelength A = 1.755 ± 0.003Aand find a K(3/Ka ratio of 0.18 ± 0.09 (90% 
confidence), consistent with the theoretical value of ~ 0.13 for low optical depth in the lines 
(Kaastra & Mewe 1993). 

The 5.36A line which we attribute to S IV- VIII is absent in either of the low-luminosity 
states, as seen in Figure 5. There is marginal evidence for a line near 5.lA in the Quiescent 
spectrum that is not present in the Ending state. This line could result from more highly 
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ionized sulfur. 

We also examine the relation between the continuum flux level and the fluorescent line 
strength. We show in Figures 6 and 7 that overall the flux in the fluorescent lines increase 
approximately linearly with the continuum flux. 

3.3. High Ionization Lines 

The low energy (higher wavelength) lines, shown in Figure 4 are seen most prominently 
outside of the Flaring period, when the low-energy continuum has been diminished. The 
lines are similar to those seen by Schulz et al. (2002b) during X-ray eclipse in the wind-fed 
X-ray binary Vela X-l. 

The flux in the low energy lines is sufficient to account for the "soft X-ray excess", a 
feature of the X-ray spectrum previously without unique explanation. HWK attributed this 
excess mainly to the passage of the X-rays through a partially ionized absorber, but they 
found that this model could not account for the excess in the spectrum at <fi > 0.6. HWK 
attributed the much larger excess seen at these phases to scattering. 

The line widths and centroid wavelengths are given in Table 4. The fluxes given are the 
observed fluxes of the lines, as they are behind an unknown column of absorbing gas. The 
errors in the line fluxes are probably overestimates, as they are based on the errors in the 
FWHM and line amplitudes, and these errors may be correlated. 

These lines do not appear to be greatly variable as a result of either changes in orbital 
phase or the intermittent flares. 

In Figure 5 we compare the low energy lines seen in the Flaring, Ending and Quiescent 
periods. We show only the spectrum during dips in the Flaring period, so that we can see 
the effect of the change in orbital phase. We have further subdivided the Ending period 
to separate periods in which the count rate was similar to the count rate in the Quiescent 
period. These time intervals are indicated by horizontal bars in Figure 1. 

The low energy lines offer several diagnostics of the emitting gas. The Helium-like 
triplets of Si XIII and Mg XI are diagnostics of density and can be used to distinguish pho- 
toionization equilibrium from collisional ionization equilibrium (Porquet et al. 2002). These 
triplets consist of a forbidden line (f, with line strength z), two closely spaced intercombi- 
nation lines (i, with strengths x and y), and a resonance line (r, with line strength w). We 
show the He-like Si triplet during the E period and the gaussian model in Figure 8. 
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The fits to the Si XIII lines (given in Table 4) give (w,x + y,z)=(10±3, 3.8^ 8 , , 5.9±™), 
where the errors are 90% confidence levels. This result is incompatible with the assumption 
of a purely photoionized plasma, which would have z > w, and instead suggests a low density 
hybrid between photoionized and collisionally ionized gas. 

Anomalies in the He-like triplets in X-ray binaries are reported in Jimenez-Garate et al. 
(2002) and Wojdowski et al. (2002). Jimenez-Garate et al. suggested from Newton-XMM 
observations of Hercules X-l that the f levels could be excited to the i levels by UV photons. 
However, this mechanism would not affect the r/f ratio. In addition, the wavelength of the 
UV photon required to raise Si XIII from the f to the i level is below the Lyman limit at 
91 1A, and would be absorbed before reaching far in the wind. 

Wojdowski et al. analyze Chandra observations of Centaurus X-3 and find similar 
results to those presented here, namely a r/f ratio near unity, inconsistent with a purely 
photoionized gas. Wojdowski et al. suggest that photoionization equilibrium still holds, but 
that resonance scattering of the continuum by the r lines adds to their flux during eclipse, 
when the direct continuum is diminished by absorption. Outside of eclipse, however, one 
expects to see the r line superimposed with an absorption feature caused by scattering. 

As can be seen in Figures 4 and 5, the ratio r/f — w/z is J> 1 even in the low part 
of the Flaring period. This is centered near = 0.65, far enough from eclipse that one 
might expect to see diminution of the r lines from an absorption trough. In §3.4, we discuss 
scattering scenarios in more detail. 

The Mg XI triplet, seen in Figure 4b, c, is in the ratio (w, x + y,z) = (6 ± 3, 3 ±2, 4 ± 2); 
unfortunately the error bars are larger than for Si, but the ratio is still consistent with a 
hybrid plasma. 

3.4. Line Formation Models 

We interpret the line strengths and ratios using the XSTAR photoionization code (Kail- 
man & Bautista 2001, Bautista & Kallman 2001). 

We used XSTAR version 2. Hi, to compute the response of a gas to incident radiation. 
The results depend most strongly on the ionization parameter £ = L/nr 2 , where L is the 
incident luminosity, n is the particle density, and r is the distance from the X-ray source, 
but also depend on the detailed source spectrum. We used a power law spectral shape with 
a high energy cut-off as determined from BeppoSax observations (Reynolds et al. 1999). 

Our line models were motivated by the suggestion, from the line ratios in He-like Si, that 
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the gas was not entirely photoionized but that dynamical heating could play an important 
role as well. Fransson & Fabian for example show that a photoionization wake can form in 
a wind-fed X-ray binary, involving a collision of ionized gas and shadowed gas accelerated 
by the normal stellar wind mechanism. Temperatures 10 6 K can form over about a stellar 
radius at an angle of 45 degrees from the X-ray source. 

Most low energy lines can be reproduced by a plasma fixed to temperature T = 5 x 10 6 K 
and with an ionization parameter log£ = 2. Such models do not predict any prominent 
RRC features, consistent with the observations. For high temperatures, RRC features are 
broadened so that they are indistinguishable from continuum (although the absence of RRC 
features does not imply high temperature.) 

However, concerns that temperatures reached in a photoionization wake might not be 
high enough to contribute to the ionization of the X-ray illuminated wind motivated us to 
study this problem further. A pure photoionization model could explain the anomalous He- 
like Si ratio through resonance scattering. Resonance scattering is a zero-sum process, and 
one might expect to have as many photons in the area scattered away as scattered toward 
the viewer. However, our point of view is special in that the inclination of the orbit is high; 
at the orbital phases of the Chandra observation, we may see the X-ray source through an 
enhanced-density gas stream. Our model of the continuum as including a Compton-scattered 
component that passes through a lower column density supports this. 

If the gas is in pure photoionization equilibrium, it would have implications for how the 
X-ray lines change with the ionizing continuum. 

Assuming pure photoionization, we have computed a Differential Emission Measure 
(DEM), a measure of the amount of gas at each bin of a parameter, usually temperature for 
work on stellar coronae in thermal equilibrium or ionization parameter £ in environments like 
gas around X-ray binaries. The problem as typically stated is overdetermined, that is, there 
are more free parameters (for example, bins of log £) than data points. However, a "penalty 
function" steers the solution toward one that is, for example, as smooth as possible. See 
Huenemoerder, Canizares, & Schulz (2001) for a more detailed description of the method. 
The major change we have introduced is to use the ionization parameter £ instead of the 
temperature T. 

When we allow values of log£ from log£ = to log£ = 5, and allow elemental abun- 
dances to be free parameters, we find that the DEM during the Ending period has two 
peaks, one at log£ < 1 (the emission of the Fe K line peaks at log£ = 1 — 2) and the other 
at log £ ~ 4. 

If the only difference between the Quiescent and Ending period is the factor of ~ 4 
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increased luminosity in the Ending period, one would expect the DEM during the Q period 
to show a peak at log£ = 3.4 instead of the log£ = 4 peak in the E period. Thus one would 
expect emission lines formed at log£ = 2 to be brighter during the Q period than the E 
period. As a result, we would expect the Si XIV line, which has peak emissivity at log£ = 3, 
to be dimmer by a factor of 2 during the Ending period than the Quiescent period, when 
in fact we have observed the line to be 2.8±0.7 times brighter in the Ending period. The 
change in the fluorescent lines between the E and Q periods is, however, consistent with the 
DEM. We find similar results when we freeze the elemental abundances to cosmic values. In 
this case, the DEM predicts that the Si XIV line would be 5 times brighter in the Q period. 

Together with the anomalous He-like Si line ratios, this might suggest the gas is not in 
pure photoionization equilibrium. However, there are large systematic uncertainties in the 
method and possible complications in the physical picture. We suggest that if the anomalous 
He-like Si line ratio results from scattering, that further observations near = 0.5 might 
reveal an absorption trough instead of emission in the Si XIII resonance line. Near = 0.5, 
the gas stream is out of the line of sight to the X-ray source so that the 6-7A continuum we 
see may be direct at not Compton scattered. The contribution of a resonance-scattered line 
component may be negligible against the direct continuum. 

3.5. Power Spectra and Quasi-Periodic Oscillations 

We created power spectra of both the overall count rate and the counts at wavelengths 
A > 4.5A, displayed in Figure 9 and Figure 10, respectively. The power spectra are all 
Leahy- normalized (Leahy et al. 1983). 

We fit the power spectrum of the integrated counts to an empirical model that consists 
of a power law with exponential cutoff, two Lorentzians, and a white noise component. 
With Leahy normalization, the white noise from counting statistics would be fixed at 2, but 
detector dead time can reduce this slightly, so the white noise level is left as a free parameter. 
We have binned the frequencies in the power spectra by a factor of 10 in order to determine 
errors in each bin. The reduced x 2 value of our model of the power spectrum is xl = 1-4- 

One Lorentzian component, detected with marginal significance, may represent a quasiperi- 
odic oscillation (QPO) with a centroid frequency of 6.5 mHz and a Q = u/Au value of ~ 15. 
This component represents a fractional rms power of 1.4%. In fact, there are several peaks 
in the 10~ 3 to 10~ 2 Hz range that may be QPO seen at low significance. Power spectral 
peaks with frequencies of exactly 1, 2 x 10~ 3 Hz probably result from the spacecraft's 1,000 
second dither period. 
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The power spectrum of the region A > 4A shows more prominent QPO features (Fig- 
ure 10). We also allow a third Lorentzian at the second harmonic of the 6.5 mHz QPO. The 
width of the second harmonic line is fixed to twice that of the fundamental. The fit gives 
xl = 2.1, and it is clear that there are bumps in the 1-10 mHz range that the model does 
not account for. The 6.5 mHz line has a fractional rms amplitude of 4.5% while its second 
harmonic has an amplitude of 4.0%. 

Similar temporal behavior has been seen from 4U 1700-37 before in the integrated 
counts observed with EXOSAT (Doll & Brinkmann 1987). Doll & Brinkmann argued that 
the timescales of variability in 4U 1700-37 reflected the flare recurrence rate and was linked 
to wind flow times. 

Examination of the cross-correlation coefficient between the hard and soft X-rays shows 
a peak consistent with zero delay between the two. 

We have also formed a power spectrum without any binning, in order to see more clearly 
the low frequency variations in the A > 4A spectrum. This power spectrum is displayed in 
Figure 11. To the eye there appears to be a cutoff near the frequency v = 3 x 10~ 4 Hz. The 
time scales of the cutoff and QPOs (10 2 — 10 4 seconds) lie between the light travel-time and 
wind-travel time between the stars. Our continuum models suggest that much of the flux in 
this wavelength range results from Compton scattering. The cutoff in the power spectrum 
could then result from smearing of the X-ray variability by Compton scattering. 

Although the region A > 4A is dominated by lines during the Ending period (Figure 5), it 
is dominated by the continuum during the Flaring period (Figure 2). Thus the low frequency 
noise and stronger QPOs in the A > 4 A region probably result from the continuum and not 
the lines. However, restricting analysis to only the wavelength regions surrounding the lines 
and to only the Quiescent and Ending periods, we still find significant low frequency noise. 

4. Conclusions and Future Work 

The current data set and the 4U 1700-37 system offer outstanding diagnostics of stellar 
winds and their disruption by compact objects. 

The low energy lines (excluding the S line at 5.36A and the Si line at 7.1lA) arise 
in highly ionized gas, yet a soft component is still visible through eclipse (HWK). XSTAR 
simulations can reproduce the low energy lines either in photoionization equilibrium with 
log£ = 2.5—3.0 or as a hybrid plasma at T = 5x 10 6 K with log£ = 1.5—2.0 (at the lower end 
of that range the Si XIVA6.18 to Si XIIIAAA6.64,6.68,6.74 ratio becomes lower than observed). 
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In Figure 12 we show that for a uniform wind model, a region in photoionization equilibrium 
with log£ = 2.5 is smaller than the primary and would not be visible in mid-eclipse. While 
this would seem to support an extended region of hot gas (perhaps a photoionization wake) 
real stellar winds are probably highly nonuniform, with density contrasts of ~ 1000. In such 
a nonuniform wind, there could be pockets of gas with log£ = 2.5 that extend all the way 
to the region indicated by log £ = for the uniform wind case. 

The Fe lines (and the low ionization SA5.36 and SiA7.11 lines) may arise in a region with 
log£ < 1. This region could be large enough for the lines to be visible through mid-eclipse 
even for a uniform wind. 

Follow-up observations can test whether there is a variable delay between the hard and 
soft X-ray oscillations. This could result, for example, if the soft X-rays arise in a more 
extended region in response to hard X-rays from the compact object. During the present 
observations, we detected no delay. 

QPO with mHz frequencies have been seen from a number of X-ray binary pulsars 
(Boroson et al. 2000). The detection of mHz QPOs in 4U 1700-37 thus suggests that the 
compact object is a neutron star. Further observations, detecting QPOs at several X-ray 
luminosity states, can can test any relationship between the frequency and an inner-disk 
radius. 

The cause of the anomalous ratio in the Si XIII triplet has not been settled, but the two 
possible explanations, resonance scattering or a hybrid plasma, would both have important 
implications. The optical depth at 6.7A, and thus density in the wind would be constrained 
by the requirement that enough photons could be scattered. On the other hand, if the line 
ratios are caused by hot gas, this contributes evidence for and helps diagnose the physical 
extent and state of a photoionization wake. 

One promising method for diagnosing resonance scattering is high resolution X-ray 
spectroscopy, but near = 0.5, as described in §3.4. 

If there is extended hot gas, we would expect it to contribute to the continuum emission. 
Such a component has been used to fit the continuum spectrum before. Haberl & Day (1992) 
used a bremsstrahlung component with kT = 0.5 keV behind a column density N# = 5 x 
10 21 cm" 2 to fit GINGA spectra (as well as the EXOSAT data of HWK). ROSAT observations 
found a bremsstrahlung component with kT = 0.47 keV and kT = 0.74 keV after eclipse 
(Haberl, Aoki, & Mavromatakis, F., 1994). Most recently, BeppoSax observations found a 
T = 0.2 ±0.1 keV component (Reynolds et al. 1999). After adjusting the normalization of 
the HEG and MEG spectra, we find that we do not require a bremsstrahlung component 
to fit the broad-band spectrum during our Chandra observations. However, our models of 
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the broad-band continuum are hampered by uncertainties in the HEG and MEG effective 
areas, and a bremsstrahlung component may be present at flux levels lower than historically 
observed. 

Other HMXB, such as Cen X-3 (Wojdowski et al. 2003) or Vela X-l (Schulz et al. 
2002) also show similar Si XIII diagnostics, but these have been interpreted to be result of 
resonance scattering. Given that Vela X-l probably also has a photoionization wake (Kaper 
et al. 1994), emission from a hybrid plasma should also be considered as a possibility, 
particularly for that system. Again, analysis of observations at a variety of orbital phases 
can help distinguish between the possibilities. Hydrodynamic models of individual systems 
can also suggest whether the wind can be heated sufficiently by dynamical means. 

SDV supported in part by NASA (NAG5-6711), and the Chandra X-ray Center (GO0- 
1100X). We would like to thank Mr. Andrew Beltz and Mr. Corey Casto for their assistance. 
We would like to thank the referee for comments. 



Table I. Parameters of the 4U 1700-37 System 



Parameter name 


Expression 


Value 


Reference 11 





Distance 


D 


1.9 kpc 


1 


Temperature of star 


T* 


42000±2000, 35000±1000 K 


2,4 


Star Radius 


Ro 


18±3, 21.9+^R 


2,4 


Mass of star 


M 


52±2,58±11 M 


2,4 


Mass of compact star 


M x 


1.8±0.4, 2.44±0.27 M© 


2,4 


Spectral Type 




06.5Iaf 


2 


Orbital inclination 


i 


> 80° 


2 


Orbital separation 


asmi 


2.0 ±0.4 R 


2 


Wind terminal velocity 


Voo 


2100±400, 1700 


2,3 


Wind mass loss rate 


M 


6, 9.5xl0~ 6 M Q yr- 1 


2,4 



a l: Ankay et al. 2001, 2: Heap k Corcoran 1992, 3: van Loon et al. 2001, 4: 
Clark et al. 2002 



Table 2. Fits to the X-ray Continuum 



Parameter 


pa 


FH 


FL 


Q 


E 


EL 


Exposure time (s) b 
N H (10 23 cm- 2 ) 
Power Law Norm c 
Power Law 7 d 
N H (2) (10 22 cm" 2 ) 
Power Law (2) Norm 


11420 

n qi+0.14 

n 1 c+0.03 
°' 16 -0.01 

1 09+° 10 
i - UM -0.07 

3.4±0.5 

o m+ 002 

1.06,607 


6164 
092 +o.i9 

u ' yz -0.12 

17+ - 04 
V- 1 ' -0.02 

1 04+ ' 13 
o fi +0.5 

06+ ' 04 
1.15,301 


5256 

1 06+° 11 
1 - uo -0.06 

n ,0+0.06 

2. Oil. 5 

002+ 006 
u.uuz_ 001 

1.13,357 


10000 
1.27±0.03 
0.038l°° 2 ° 

0.29, < 3.9 
1.0 ±0.4 x 10~ 4 
1.09,61 


21531 
2.03±0.08 
0.095±0.018 

l.ooio.io 

0.95 (< 1.8) 
3l; 2 x 10~ 3 
1.07,301 


3300 

1 8+0-3 

04+ ' 03 
u ' u4 -0.02 

1 2+ - 4 
0.97 (< 2.1) 
4±2 x 10" 3 
(c-stat?61) 



a F=Flaring, FH=Flaring High, FL=Flaring Low, Q=Quiescent, E=Ending, EL=Ending Low, as shown in Fig- 
ure 1 

b The total on-source time, including non-readout time ("dead time") 
c The normalization is denned as the photons keV -1 cm~ 2 s _1 at 1 keV. 

d This is the "photon index", so that the number of photons at energy E is given by K(E/lkeV) -7 , where K is 
the normalization 

c For the low ending state, we used the c-statistic (Cash 1979) instead of x 2 , as even when the wavelength bins 
were grouped, the statistics were too poor to apply x 2 . 
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Table 3. Measurements of Line Wavelengths, Widths, and Fluxes 



A (A) 


Mode a 


Ion 


A ( , (A) b 


FWHM (km/s) 


Flux (10 6 photons s 1 cm 2 ) 




1.753±0.001 


F 


FeII-XII,K/3 c 


1.756 


2200±700 


300±200 


1.76±0.01 


FH 






2000±900 


300±300 


1.75±0.01 


FL 






2000±900 


200±200 


1.755±0.003 


E 






900±400 


200±100 


1.937±0.002 


F 


FeII-XII,Ka 


1.936-1.937 


2200±700 


1200±200 


1.937±0.002 


FH 






2000±900 


1400±300 


1.937±0.002 


FL 






2000±900 


900±200 


1.940±0.002 


Q 






< 2200 


170±60 


1.936±0.001 


E 






900±400 


1100±90 


1.938±0.001 


EL 






unresolved 


350±130 


3.355±0.001 


E 


CaVI-VII 


3.356-3.354 


unresolved 


46±15 


4.18 d 


F 


Ar VI-IX 


4.186-4.178 


unresolved 


40 ±30 


4.18 d 


FL 






unresolved 


66±40 


4.18 d 


Q 






unresolved 


< 13 


4.18 d 


E 


Ar VI-IX 


4.186-4.178 


unresolved 


20±9 


4.18 d 


EL 






unresolved 


6IJ -20 


4.744±0.004 


F 


SXVI 


4.7274-4.7328 


unresolved 


60±30 


4.744±0.002 


FH 






unresolved 


90±50 


4 69+ ' 05 


FL 






unresolved 


30±30 


4.74 d 


E 






unresolved 


14 ±7 


4.74 d 


EL 






unresolved 


< 22 


5.10±0.01 


Q 


SXIII-XV 


5.075-5.161 


undetermined 


U±? 


5.36 d 


F 


S IV- VIII 


5.370-5.356 


< 3000 


50±30 


5.366±0.005 


FH 






< 4000 


80±50 


5.36 d 


FL 






unresolved 


40+ 40 
4U -30 


5.37±0.03 


Q 






> 1000 


12±12 


5.359±0.005 


E 






3000±1000 


70±15 


5.36±0.02 


EL 






3000 




6.17 


F 


Si XIV 


6.199 


unresolved 


29±10 


6.17 


FH 






unresolved 


30±16 


6.22±0.01 


F 


? 




10000 


30±14 


6.22±0.01 


FH 


? 




10000 


47±25 


6.17 


FL 






unresolved 


27±11 


6.18±0.01 


Q 






1800±£°° 


6+ 4 

D -3 


6.177±0.003 


E 






1200±300 


17±4 


6.179±0.001 


EL 






1200 


13±\° 


6.62±0.01 


F 


Si XIII 


6.647 


1900±500 


38±12 


6.61±0.02 


FH 






1200±500 


39±19 


6.633±0.006 


FL 






1000±500 


25±10 


6.63±0.01 


Q 






2000±1000 


H±! 


6.640±0.005 


E 






1200 ± 400 


10 ±3 


6.640 d 


EL 






1200 




6-68±g;g5 


Q 


Si XIII 


6.684,6.687 


1400+ 1500 


1.91- 


6.684±0.009 


E 






1200±400 c 




6.684 d 


EL 






1200 


< 3 


6.73±0.01 


F 


Si XIII 


6.739 


1900±500 


33±12 


6.73±0.01 


FH 






1200±500 


43 ±20 


6.73±0.01 


FL 






1000±500 


11+? 


6.73±0.02 


Q 






2500±«S° 


6+ 4 

D -3 
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Table 3 — Continued 



A (A) 


Mode a 


Ion 


Ao (A) b 


FWHM (km/s) 


Flux (10~ 6 photons s" 1 cm" 2 ) 


6.733±0.007 


E 


Si XIII 


6.739 


1200±200 


r- n + 2 6 


6.733 d 


EL 






1200 


< 8 


7.05±0.06 


F 






1 onon+lOOOO 

13000I 4000 


50±20 


7.00±0.11 


FH 


Si blend 






120±40 


7.09±0.04 


FL 






8000I 3000 


11 ±15 


7.02±0.03 


E 


Si blend 




8000tf ° ° 


12 ±4 


7.099±0.004 


F 


<SiVII 


> 7.063 


1100 ± 400 


40±10 


7.099±0.007 


FH 






1100+ 700 


52±16 


7.098±0.002 


FL 






800±400 


32±10 


7.110±0.007 


Q 






1700™ 


14±5 


7.110±0.003 


E 






1300±300 


22 ±4 


7.110±0.01 


EL 






1300 


13±7 


7.84±0.02 


E 


MgXI 


7.850 


2100+ 1800 
ziuu_ 70() 


3±2 


8.41 f 


FH 
FL 

Q 
E 


MgXII 


8.419,8.424 


< 1400 
1100±500 

< 8000 
1300±700 


14±7 
13±6 
6±3 
5±2 


9.1681 


FH 
E 


MgXI 




3000+ 3000 
csuuu_ 1000 

2500+ 2 """ 


21±10 
4±2 


9.2280 


FH 
E 


MgXI 




o 000 +3000 
2500+ 2 """ 


< 10 
3±2 


9.3134 


FH 
E 


MgXI 




3000±?g88 
2500i 2 8 r 


< 14 
6±3 


10.238 


FH 
E 


NeXLy/3 




unresolved 
unresolved 


7 +8 
'-5 

3±2 


12.132 


FH 

E 


NcXLya 




unresolved 
unresolved 


O+10 

4±3 



a F=Flaring, FH=Flaring High, FL=Flaring Low, Q=Quiescent, E=Ending, EL=Ending Low, as shown 
in Figure 1 

b Rest wavelength. 

c We give wavelengths for Kai. The centroids are « 0.002A greater than the a\ wavelengths. The fluxes 
we give are for Kai+Ka2 = 1.5 K Q1 

d Marginal detection. The wavelength was frozen to the given value 

°The He-like Si lines are fixed to have equal FWHM. 

f All the lines with A > 8A have their wavelengths frozen to the given value 
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Fig. 1.— The combined MEG and HEG light curve of 4U 1700-37. We have divided the 
light curve into three separate intervals, which we call Flaring, Quiescent, and the Ending 
periods. To separate the effects of orbital phase and X-ray flaring, we investigate times in 
the F period when the X-ray light curve, smoothed over 30 pixels, is less than 30 cts per 
3.24 second readout bin. We also investigate the spectrum during times in the E period in 
which the count rate resembles the count rate during the Q period. The separate intervals 
in the F and E periods that are selected for further analysis are indicated by the horizontal 
lines with bars at the ends. We overplot in red the contemporaneous RXTE light curve. The 
scale of the RXTE light curve is given on the right y axis. 

Fig. 2. — The combined MEG and HEG spectrum of 4U 1700-37 integrated over the Flaring 
period (F), Quiescent period (Q), and Ending period (E) of our observation. The 1.93A Fe K 
line is apparent, as well as lines in the 4-8A range in the F and E spectra. The spectra have 
been binned to 0.04A resolution. We also show our model of the continuum, which consists 
of a power law behind two different column densities. The softer, less-absorbed spectrum 
is plotted with dashed lines, while the harder, more-absorbed component is plotted with 
dotted lines. The colors are chosen consistently within each observing segment (F,Q,E). 

Fig. 3. — The Fe Ka and K/3 lines during the (a) Flaring (0.005 A/pixel), (b) Quiescent 
(0.01 A/pixel), and (c) Ending (0.005 A/pixel binning) periods. The solid curve is the model 
fit, while the histogram is the combined MEG and HEG data. The edge near 1.74A results 
from neutral iron and is implicit in the model for cold absorption that the continuum (power 
law) passes through. 

Fig. 4. — The 1.5-13A spectrum during the Ending period, with Gaussian fits to the lines 
and line identifications. We show (a) the count spectrum, (b) the log of the count spectrum, 
(c) the spectrum with the continuum model subtracted. 

Fig. 5. — A comparison of the MEG spectra (in counts per second per 0.02 A bin) at 
different times. From top to bottom, the panels show: the spectrum during the low count 
rate intervals of the Flaring period, the spectrum integrated over the entire Quiescent period, 
the spectrum integrated over the entire Ending period, and the spectrum during the low 
count rate intervals of the Ending period. The error bars show the la errors for the 7.1 A 
line peak. The times of the low count rate periods are indicated by the horizontal bars in 
Figure 1. 

Fig. 6. — Correlations between the flux in the X-ray continuum (1-13AA) and the flux mea- 
sured in fluorescent lines. For the continuum flux, we include both the direct and Compton- 
scattered components at the periods (FH,FL,Q,E,EL). a),b),c),d), and e) denote whether 
the fluorescence feature is due to Fe K, Ar, S, narrow Si near 7A, or broad Si near 7A. All 
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error bars are 90% 

Fig. 7. — This is identical to the previous plot, but now the continuum includes only the 
hard component but without any absorption. Again, the fluorescence feature is caused by 
a) Fe K, b) Ar, c) S, d) narrow Si near 7A, or e) broad Si near 7A. 

Fig. 8. — The spectrum near the He- like Si-triplet in the Ending period. The rest wavelengths 
of the Resonance, two Intercombination lines, and the Forbidden lines are shown with vertical 
dotted lines. 

Fig. 9. — The Leahy- normalized power spectrum of the continuum. We have subtracted 
counts in the low energy lines from the total counts in HEG and MEG 1st order. The total 
model and individual model components are also shown. We show error bars where the 
logarithmic scale allows them to be visible. 

Fig. 10. — The Leahy-normalized power spectrum of the low energy region. The total model 
and individual model components are also shown. 

Fig. 11. — The Leahy-normalized power spectrum. In order to emphasize the low-frequency 
behavior, we have not performed any binning. 

Fig. 12. — Contours of constant ionization parameter log£ for a simple model of a uniform 
wind with M = 6 x 10~ 6 M yr _1 , a terminal velocity of = 1700 km/s, a wind velocity 
law of v(R) = i>oo(l — Ro/R) 13 where (3 = 0.8. The "x" marks the position of the compact 
object and the circle shows the primary. The log£ = 0.0 ionization boundary is cut off where 
the primary's shadow prevents X-rays from reaching the wind. 
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